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["T , I ABSTRACT 

We consider mean- field dynamo models with fluctuating a effect, both with and with- 
out large-scale shear. The a effect is chosen to be Gaussian white noise with zero 
mean and a given covariance. In the presence of shear, we show analytically that (in 
infinitely large domains) the mean-squared magnetic field shows exponential growth. 
The growth rate of the fastest growing mode is proportional to the shear rate. This 
result agrees with ea rlier numerical res u lts of lYousef et al.l (|2008|) and the recent an- 
alytical treatment bv lHeinemann et al.l ([201 ll ) who use a method different from ours. 
In the absence of shear, an incoherent a 2 dynamo may also be possible. We further 
show by explicit calculation of the growth rate of third and fourth order moments of 
the magnetic field that the probability density function of the mean magnetic field 
generated by this dynamo is non-Gaussian. 
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The dynamo mechanism that generates large-scale magnetic 
fields in astrophysical objects is a topic of active research. 
Almost all astrophysical bodies, e.g., the Galaxy, or the 
Sun, show presence of large-scale shear (differential rota- 
tion). It is now well established that this shear is an es- 
sential ingredient in the dynamo mechanism. This view is 
also supported by direct numerical simulations (DNS). In 
particular, DNS of convective flows have been able to gen- 
erate l arge-scale dynamos predominantly in the pres ence of 
shear JKapyla et all 120081 ; iHughes fc Proctoil l2009l h while 
non-shearing large-sc ale dynamos are on ly possible at very 
high rotation rates l)Kapvla et al.l I2009T ). The other vital 
constituent of the large-scale dynamo mechanism is helic- 
ity of the flow which is often described by the a eff ect. In- 
deed most dynamos, including the early model by IParkerl 
(1955), are the result of an a effect combined with shear. 
Of these two ingredients, shear is typically constant over 
the time scale of generation of the dynamo. For the case 
of the solar dynamo, shear or differential rotation are con- 
strained by helioseismology. By contrast, measuring the a 
effect is a non-trivial exercise. Whenever it has been ob- 
tained from DNS studies, it was f ound to have large fluctua - 
tions in space and time; see e.g.. IBrandenburg et al.l (2008) 
for the PDF of diff e rent c omponents of the a tensor and 
ICattaneo fe Hughea (|199a ) for the fluctuating time series 
of the total electromotive force, which is however different 
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from the a effect. What effect do these fluctuations have on 
the properties of the a-shear dynamo? In particular, can a 
fluctuating a effect about a zero mean drive a large-scale 
dynamo in conjunction with shear? Recent DNS studies 
JBrandenburg||2005al : lYousef et af]|2008t iBrandenburg et all 
I200SI . hereafter referred to a s BRRK) sug g est th at the an- 
swer to this question is yes. lYousef et al.l (20081 ) have fur- 
ther provided compelling evidence that the growth rate of 
the large-scale magnetic energy in such dynamos scales lin- 
early with the shear rate, and the wavenumber of the fastest 
growing mode scales as the square root of the shear rate. 
It has been proposed that such a dynamo can also emerge 
by an alternate mechanism (which has nothing to do with 
fluctuations of the a effect) involving the interaction be- 
tween shear and mean current density. This mechanism 
is ca lled the shear-current effect (|Rogachevskii fc Kleeorinl 
I2003T) . H owever recent DNS st udies of BRRK as well as 
those of IBrandenburg! (|2005bl ) have not found evidence 
in support of it, and analytical works (JRadler fc Stepanov 
20061 : iRiidiger fc Kitchatinovll2006l : ISridhar fc Subramania" 



2009), for small magnetic Reynolds number, have doubted 
its existence. Furthermore, the shear-current effect does not 
produce the observed scaling, namely, growth-rate scales 
with the shear rate squared, and the wavenumber of the 
fastest growing mode scales linearly with the shear rate; see 
Section 4.2 of BRRK. 

It behooves us then to consider the interaction be- 
tween fluctuating a effect a n d she ar as a possible dy- 
namo mechanism. iKraichnanl l| 19761) was first to propose 
that fluctuations of kinetic helicity can give rise to nega- 
tive turbulent diffusivities, thereby giving rise to a dynamo 
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that is effective only in small scales (|Moffattlll97a ). How- 
ever, more relevant in the pre sent context is the work of 
IVishniac fc Brandenburg! (|l997f ). who have investigated (us- 
ing numerical tools) an one-dimensional model and (with 
analytical techniques) a simple mean-field model in a one- 
mode truncation with shear and fluctuating a effect. These 
models they referred to as incoherent a-il dynamos. By di- 
viding the poloidal field component by the toroidal one, they 
turned the mean-field dynamo equations with multiplicative 
noise into one with additive noise of Langevin type. They 
drew an analogy between the behaviour of the mean mag- 
netic field and Brownian motion in the sense that the mean 
field does not grow although the mean-square field grows 
with a growth rate that scales with shear rate to the 2/3 
power. 

Over the last decade, the problem has be en stud- 
ied u sing a variety of different appro a ches (ISokoloyl 
19971; ISilanFevI l20od; iFedotoy et all l200fj; IProctorl l2007t 



Kleeorin fc Rogachevskiil |200S| ; ISur fc Subramanianl I2009D . 
Models with fluctuating a effect, which depend on both 
space and time, can be divi ded into two ca t egories, one in 
which q is inhomogeneous (ISilant'evi |2000| ; iProctorl 120071 ; 
iKleeorin fc Rogachevskii 2008) and the other in which it is 
a constant in space. (There are other technical differences in 
the two approaches as well.) In this paper we shall confine 
ourselves to the second class of m odels. Fo r the fi rst cate- 
gory, using a multiscale expansion, IProctorl (|2007l ) found a 
quadratic dependence of growth rate on shear rate. Th is re- 
sult is contradicted by Kleeo rin fc Rogachevskiil (J2008J) who 
found that the mean field averaged over a does not show a 
dynamo and that the fluctuating a effect adds to the diag- 
onal components of the turbulent magnetic diffusivity ten- 
sor, in which they consider the first and second moments 
(calculated over the distribution of a) of the horizontally 
averaged mean field. The first moment shows the absence 
of dynamo effect, but the second moment shows exponen- 
tial growth with t he same scaling behaviour observed by 
lYousef et al.l ([2008). Motivated by their study, we solve here 
a similar model using a technique known as Gaussian in- 
tegration by parts to obtain the same results. We further 
show that, even in the absence of shear, it may be possi- 
ble for incoherent a 2 dyn a mos t o operate. For the model of 
IVishniac fc Bran denburg! (|1997l ) we can also calculate the 
higher (third and fourth) order moments of the magnetic 
field to demonstrate that the probability distribution func- 
tion (PDF) of the mean magnetic field generated by an in- 
coherent a-shear dynamo is non-Gaussian. 



2 MEAN-FIELD MODEL 

O ur mean-field model is designed to describe the simulations 
of lYousef et al.l (|2008r ) and BRRK. In particular, there is a 
large-scale velocity given by U = Sxe y . The turbulence is 
generated by an isotropic external random force with zero 
net helicity. The mean fields are defined by averaging over 
two coordinate directions x and y. By assuming shearing- 
periodic boundary conditions, as done in the simulations, 
this averaging obeys the Reynolds rules. By virtue of the 
divergence- less property of the mean magnetic field, v.i.z., 
djBj = 0, and in the absence of an imposed mean field, B z 



can be set to zero. The resultant mean-field equations then 
have the following form (BRRK) 

d t B X — —OLy X d Z B X — CtyydzBy — T)y X d z By + T]yyd z B X , (l) 

d t By — SB x +a xx d z B x +a xy d z By~ri X yd z B x +rj xx d z By. (2) 

Here, B x and B y are the x and y components of the mean 
(averaged over x and y coordinate directions) magnetic field, 
ntj are the four relevant components of the turbulent mag- 
netic diffusivity tensor and ohj are the four relevant compo- 
nents of the a tensor. The shear-current effect works via a 
non-zero r) yx , provided its sign is the same as that of S. In 
this paper we do not consider the possibility of the shear- 
current effect; hence we set rj yx — 0. Here we have ignored 
the molecular diffusivity, as we are interested in the limit of 
very high magnetic Reynolds numbers. 

The fluctuating a effect is modelled by choosing each 
component of the a tensor to be an independent Gaussian 
random number with zero mean and the following covariance 
(no summation over repeated indices is assumed): 

(o«(t)a«(0> = D ij 5 ik 5 j i5{t - t'). (3) 

We also assume that the ctijS are constants in space. Note 
that DNS studies of BRRK have shown that the coefficients 
of turbulent diffusivity also show fluctuations in time, but we 
have ignored that in this paper. To make our notation fully 
transparent, let us clearly distinguish between two differ- 
ent kinds of averaging we need to perform. The mean fields 
themselves are constructed by Reynolds averaging, which in 
our case is horizontal averaging, and is denoted by an over- 
bar. As we are dealing with mean field models, the quanti- 
ties appearing in our equations have already been Reynolds 
averaged. But, as the ajj in our mean field equation are 
stochastic, we study the properties of our model by writ- 
ing down evolution equations for different moments of the 
mean magnetic field averaged over the statistics of <Xij ■ This 
is denoted by the symbol (•). To give an example, the first 
moment (mean over statistics of a) of the x component of 
the mean magnetic field is denoted by (B x ) and the second 
moment (mean square) is denoted by {B x }. 

In numerical studies, such averaging has to be per- 
formed by averaging over sufficiently many simulations 
with independent r e alisatio ns of noise, as was done by 
ISur fc Subramanianl (J2009J). In controlled experimental 
setup such ensemble averaging is done by executing the ex- 
periment several times with idential initial conditions. In na- 
ture we expect the noise to have a finite correlation length. 
Averaging over different realisation of the noise can then be 
done by averaging over spatial domains of size larger than 
the correlation length of the noise. Note further that in the 
present case the components of B undergo a random walk 
and are hence non-stationary. Therefore we cannot replace 
ensemble averages by time averages. 



3 RESULTS 

3.1 First and second moments 

It is convenient to study the possibility of a dynamo effect 
in Fourier space. Under Fourier transform, 



B x 



B x e- lkz dk, By 



(4) 
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Equations |T} and ([2]) transform to, 

Ot-Dx = Oiy X %KD X OLyylKlDy + T)y X K Dy ^] yy K ^X 



(5) 



dtBj, = SB X + a xx ikB x +a xy ikB y +n xy k B x ~rj xx k B y . (6) 

Equations ((5]l and (J6]) are a set of coupled stochastic differ- 
ential equations (SDEs) with multiplicative noise. 

We now want to write down an evolution equation for 
C 1 = {{B x }, (By)), which is the mean field averaged over the 
statistics of otij. Note that C 1 is the first order moment of 
the probability distribution function (PDF) of the magnetic 
field. As the a effect is taken to be Gaussian and white-in- 
time it is possible to obtain closed equations of the form 

dtC 1 = NiC 1 

with 

\Hyy ' y x ) ^y x 

o T" k T\ X y k \f] xx -\- D X y) 



Ni 



(7) 



(8) 



If there is no shear-current effect, i.e. r/ yx — 0, there is no dy- 
namo. The fluctuations of the a effect actually enhance the 
diagonal components of the turbulent magnetic diffusivity. 
This result agrees with iHeinemann et al.l l)201ir ) who used 
a different mo del and found that the mean ma gnetic field 
does not grow. IVishniac fc Brandenburg! (119971 ) also found 
the same for their simplified zero-dimensional model. How- 
ever, if one assumes that different components of the a ten- 
sor are correlated, a different result can be obtained; see 
Section 13.31 

We have used two different technique s to derive the ma- 
trix N i in ((SJ). In Appendix [XI following iBrissaud fc Frischl 
(|1974T ). we have used a perturbation expansion in powers 
of noise strength. This expansion works even when a is not 
white in time. However, for white-in-time a it is enough to 
retain only the leading order term. In Appendix [B] we have 
used Gaussian integration by parts which works because of 
the white-in-time nature of the a effect. This method can 
be easily applied to study higher order moments of the PDF 
of the magnetic field too. Hence we shall use it extensively 
in the rest of this paper. 

Although the mean magnetic field does not grow, the 
mean-squared magnetic field can still show growth. To study 
this we now write a set of equations for the time evolution of 
the second moment (covariance) of the mean magnetic field 
averaged over ohj . We emphasise that we are not considering 
here the covariance of the actual magnetic field that would 
include also the small-scale magnetic fluctuations, whi ch are 
relevant to the small-scale dynamo (|Kazantsevl ll968). The 
covariance of the actual magnetic field was later also con- 
sidered by Hovng ( 1987) in c onnec tion with a-Q. dynamos. 
Following IHeinemann et al.l ( J2011T ) , we define a covariance 
vector, 



C 2 = ((B X B* X ), {B y B' v ), {B x B* y ) + (B:6 v )). 



(9) 



The evolution equation for the covariance vector is given by 



d t C 2 = N 2 C 2 , 
where N2 is given by 

2k 2 D y 



(10) 



2k T\yy 

f)xy) 



2{S + k 2 



£™ r/ xx 



Vy x 
D -f- K Tj X y 
-K \Lfy x ~Y 1-J x y I ^rj xx J 



The only non-trivial terms in the derivation of the above 
equation are the terms which are products of components 
of q effect and two components of the magnetic field. We 
evaluate them by using the same technique used to obtain 
(0 and ((HJ; see Appendix [B] for details. 

The characteristic equation of the matrix N2 is a third 
order equation, the solutions of which gives the three so- 
lutions for the growth rate 27. For simplicity let us also 
choose n xx = rj yy = rj t . In other words, we take the turbu- 
lent magnetic diffusivity tensor to be diagonal and isotropic. 
We further note that D yx and D xy contribute only in en- 
hancing the turbulent magnetic diffusivity of C3, so we can 
safely ignore them compared to rj xx . With these simplifying 
assumptions the equation for the growth rate reduces to 



C 3 - 4:k 2 D yy [S 2 + k 2 D xx (\ = 0, 



(11) 



where £ = 2(k 2 nt + 7). For large enough S we can always 
ignore the second term inside the parenthesis of (|lip . This 
gives the three roots of 7 as 



-k\ + [\k 2 DyyS 



, !/3 

(l,w,w ) 



(12) 



where (1, w, uj 2 ) are the three cube roots of unity, of which uj 
and u> 2 have negat ive real parts. The same dispersion rela- 
tion is obtained by IHeinemann et al.l |201lr ). The wavenum- 
ber of the fastest growing mode, fc pcak , is given by 



k 



peak 



,1/2 



(Ei 



\^4 



1/4 



(13) 



The growth rate of the fastest growing mode is given by 

1/2 



"£('-£)fr 



ISI. 



(14) 



This is th e sam e scaling numerically obtained by 
lYousef et ail (|200Sl ). 



3.2 Incoherent a dynamo 

Let us now consider a different case where shear is zero. In 
that case, (|11[) becomes a quadratic equation in £, 



e 



4k*D XX Dyy = 0, 



with solutions, 



(-1* ± v 7 ^ 



:J-Sy 



(15) 



(16) 



Hence, it may be possible for fluctuations of a to drive a 
large-scale dynamo (in the mean-square sense) even in the 
absence of velocity shear. 

To summarise there are two possible dynamo mecha- 
nisms in our dynamo model. In both of them the magnetic 
field grows in the mean-square sense. The first one is an in- 
coherent a-shear dynamo. For large enough shear this is the 
fastest growing mode. However, this dynamo has no oscillat- 
ing modes because the modes for which 7 have a non-zero 
imaginary part have negative real part. An incoherent a 2 
dynamo mechanism also exists in this model. The condition 
for excitation of a fluctuating a-shear dynamo is 



(k 2 D yy S 2 /2)^ 
k 2 n t 



>1, 



(17) 
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Figure 1. Sketch of the dispersion relation, JTT} for differ- 
ent values of velocity shear S (continuous lines with difference 
colours/grey shades from bottom to top S = 0.5, 1., 1.5, 2., 2.5), 
and l|16jl (broken line) at the very bottom. Velocity shear and 7 
have dimensions of inverse time and are measured in the units of 

^rmsn-f - 



and the condition for excitation of a fluctuating a dynamo 

is 



Vt 



> 1. 



(18) 



The condition that a fluctuating a dynamo is preferred 
compared to an a-shear one is 



S 4 



>1, 



(19) 



or V2k 2 D xx /\S\ > 1 for D xx = D yy . 

To compare with DNS we need to use some estimates 
of r/t, D xx and D y v . We use r/ t = u Iuls /3kf, as obtained 



by ISur et al.l l|2008l ) without shear. A slightly larger value 
was found by BRRK in the presence of shear. We further 
use D xx — Dyy = Mrma/9. For this choice of parameters 
the incoherent a 2 dynamo does not grow. Here, u rms is the 
mean-squared velocity and fcf corresponds to the character- 
istic Fourier mode of the forcing if the turbule nce has been 
maint ained by an external force, as done by lYousef et al.l 
(2008) or BRRK. For turbulence maintained by convection, 
k{ should be replaced by the Fourier mode corresponding 
to the integral scale of the turbulence. Typically, mean-field 
theory applies for modes with k < k{. Lengths are measured 
in units of 1/fcf and velocity is measured in the unit of u rms . 
This makes l/u lrns kf the unit of time. The two dispersion 
relations are plotted in Fig. [T]for different values of S. 



3.3 Effects of mutual correlations between 
components of the a tensor 

So far we have assumed that only the self-correlations of the 
components of the a tensor are non-zero and the mutual 
correlations zero. Let us now generalise ([3]) to 



<a«(t)a H (t')>=*>S*(t-0- 



(20) 



Obviously, TT^ = T>fj and T)V. — D t j. It is again possible to 
write a closed equation for the first moment of the magnetic 
field in the form dtC 1 = NiC 1 with 



Ni = 



where 



</,r 



= -v 



k 2 {vt+Vyy) 

O T" rZ TJxy 


rC T]yx 

-k 2 (r} t + fj xx ) 


i 


yx j~*yy 
yx ^xxi 

xx > ^ yx i 


r) — -V yx +V yy 

'fyx — u yy \ ^xyi 

fi _ fjxy _ T>xx 
ifxx — ^ xy J ^yy • 



(21) 



(22) 
(23) 



This is a generalization of ([8]). Here, for simplicity, we have 
assumed r\ xx = r\ yy — r\ x and rj yx — r\ xy = 0; in other words 
the conventional shear-current effect is taken to be zero. 

Note that, unlike the self-correlation terms, i.e. T>l J -, 
which must be positive, the mutual correlation terms (e.g., 
T>xx> or T^yy) can have either sign. Hence, it is a-priori not 
clear from (|21[) whether a dynamo is possible or not. How- 
ever, we note two interesting possibilities below. First, the 
mutual correlations of the fluctuating a now contribute to 
off-diagonal components of the turbulent magnetic diffusiv- 
ity tensor. A dynamo is possible if 

JV]yx 



k 2 r/ t 



> 1. 



(24) 



Such a dynamo might look deceptively similar to the shear- 
current dynamo, but is actually not so because in the regular 
shear-current effect r\ yx emerges due to the presence of shear 
and hence must be proportional to S for small S. Hence 
for a regular shear-current dynamo we would have 7 ~ S 2 
and fc pcak ~ S; see BRRK. But here fj xy emerges due to 
fluctuations of a and may be independent of S, at least for 
small S. This would imply that 7 ~ S and fc poak ~ VS but 
this time even for the first moment of the magnetic field. 
The growth rate of such a dynamo may however be quite 
small as it is proportional to r\ yx which is the difference 
between two terms each of which are correlations between 
different components of the a tensor. Secondly, HT> XX > T> yx 
the fluctuating a effect gives negative contributions to even 
the diagonal components of turbulent diffus ivity, which is 
reminiscent of the result of lKraichnanl II 19761 ) . 



3.4 Scaling in a simpler zero-dimensional model 

The essential physics of (Jl| and ([2} can be captured 
by an even simpler mean-field model in a one-mode 
truncation, but with fluctuati ng a effect, introduced by 
IVishniac fc Brandenburd II 19971 ). Their model, rewritten in 
our notation and setting all k factors to unity is 



d t B x 
d t B v 



aB v 



SB, 



rjtB x , 

-VtBy 



(25) 
(26) 



This model can be analysed in exactly similar ways. By con- 
struction, this model does not have a fluctuating a 2 effect, 
and a is a Gaussian random variable with zero mean and 



covariancc 



{a(t)a(t')) =D5(t-t'). 



(27) 



For this model we adopt a more general framework and de- 
fine the growth-rate of the p-th order moment of the mag- 
netic field (the first order is the mean and the second order 
is the covariance ) to be P7 P . Explicit calculations, shown in 
Appendix [C] give 



7i 



-T)t, 72 



-Vt + 



Uds 2 



1/3 



,9 



2/3 



(28) 
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Figure 2. £ p versus p as obtained from (1306 ■ If the magnetic 
field had obeyed Gaussian statistics, f p versus p would have been 
constant. 



We show in Appendix [C] that this two-third scaling with 
shear rate in this zero-dimensional model is equivalent to 
7 ~ S scaling for JT} and © . 



3.5 Possibility of intermittency 

We note that (0 and (0 can be considered as coupled 
stochastic differential equations of the Langevin type but 
with multiplicative noise. We have taken the probability 
distribution function (PDF) of the noise to be Gaussian. 
But, by virtue of multiplicative noise, the PDF of the mag- 
netic field may be non-Gaussian. We have already calculated 
the first and second moments of this PDF. To probe non- 
Gaussianity we need to calculate the higher order moments. 
For ((5} and (|6} this is a formidable problem. But i t is fa r 
simpler f or the model of IVishniac fc Brandenburg |l997|) . 
ISokolovl i|1997l ) has already argu ed that the statistics of 
the m agnetic field in the model of Vishniac fc Brandenburg 
(1997) is intermittent; see also ISur fc Subramanianl (2009). 
In Appendix [C] we show that the growth rate for the third 
and fourth order moments are given by 



73 



-Vt + 



( 18DS 2 
\ 27 



1/3 



74 



-Vt + 



( MPS 2 
I 64 



1/3 



(29) 



Clearly, 7 P has the same scaling dependence on D and S, 
independent of p, but nevertheless they are different, i.e., 
the PDF is non-Gaussian. This non-Gaussianity is best de- 
scribed by plotting 



Cp = 



7p+ ?7t 

(DS- 



2-11/3 



(30) 



versus p in Fig. [5] 

Let us now conjecture that as p — > oo, £ p remains finite. 
Remembering that £i = 0, the general form would then be 



G 



(p- l)(a + axp + a 2 p 2 



1/3 



(31) 



Substituting the form back in (I28|l and (|29|l we find ao = 36, 
ai = —30, and a 2 = 7. This formula is also plotted in Fig. [2] 



4 CONCLUSIONS 

In this paper we have analytically solved a mean-field dy- 
namo model with fluctuating a effect to find self-excited 
solutions. We have studied the growth rate of different mo- 
ments (calculated over the statistics of a) of the magnetic 
field. There are three crucial a spects in which our results, 
th e DNS of Yousef et al.l (I200&I) . and the analytical results 
bv lHeinemann et al.l ( 20111 ) agree: (a) There is no dynamo 
for the first moment of the magnetic field, (b) the second 
moment (mean-square) of the magnetic field shows dynamo 
action, and (c) the fastest growing mode has a growth rate 
7 ~ 5 at Fourier mode fc pcak ~ \^S. We have further 
shown that these aspects of our results can even be repro- 
duced by a simpler zero-dimensi onal mean field model due to 
IVishniac fc Brandenburg (|1997l ). For this simpler model we 
have also calculated the growth rate for third and fourth or- 
der moments and we have explicitly demonstrated the non- 
Gaussian nature of the PDF of the magnetic field. Given 
that the incoherent a-shear dynamo (often with an addi- 
tional coherent part) is the most common dynamo mech- 
anism our results provide a qualitative reasoning of why 
large-scale magnetic fields in the universe may be intermit- 
tent. However note also that we have merely shown that 
the growth rates of the different moments of the magnetic 
field are different. The eventual nature of the PDF of the 
magnetic field will also be influenced by the saturation of 
this dynamo which is outside the realm of this paper. We 
have also shown that it is possible to find growth of the first 
moment (mean) of the magnetic field if mutual correlations 
between different components of the a tensor are assumed to 
have a certain form (Section 13. 3[1 , It will thus be important 
to check such assumptions from future DNS. 

As our paper has been inspired by iHeinemann et al.l 

(J2011T ) it is appropriate that we compare and contrast our 
model and techniques with theirs. Their model consists of 
the equations of magnetohydrodynamics (MHD) with an ex- 
ternal Gaussian, white-in-time force (in the evolution equa- 
tion for velocity) with the additional assumption that the 
non-linear term in the velocity equation is omitted. The 
model thus applies in the limit of Reynolds number Re <C 1. 
They perform averaging over a;j/-coordinates to obtain mean 
field equations with an a effect which depends on the he- 
licity averaged over coordinate directions. Our mean field 
model is derived by first averaging over coordinate directions 
(standard Reynolds averaging) with the additional assump- 
tions on the statistics (Gaussian, white-in-time) of a. Our 
results are thus not limited by the smallness of the Reynolds 
number, although all the usual limitations of mean-field 
theory apply. The assumption of the Gaussian nature of 
a is we ll supported by numeric al evidence; see Fig. 10 of 
BRRK. IHeinemann et al.l (|201ll ) have further used a quasi- 
two-dimensional velocity field, but this we feel is not an im- 
portant limitation. They average the first and second mo- 
ment of the magnetic field over the realisations of force by 
using cumulant expansion in powers of the Kubo number. 
As they truncate the expansion at the lowest order in Kubo 
number it applies to the case of small Kubo number. The 
most restrictive assumption in our model is the assumption 
of the white-in-time nature of the a effect. This assumption 
however allows us to obtain closed equati ons for all the mo- 
ments of the magnetic field. The results of lHeinemann et alj 
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i|201ll ) is not limited by this assumption. It is interesting 
to note that, even under the assumption of the white-in- 
time natur e of the a effect, we obt ain the same scaling be- 
haviour as Heinemann et al.l l|2011f ) and the DNS studies of 
lYousef et all l|2008l) . 

Here, let us also mention that 

iKolokolov. Lebedev and Sizovl (|2011T ) have recently ap- 
plied similar techniques to study small-scale kinematic 
dynamos in a smooth delta-correlated velocity field in the 
presence of shear to find 



T™=^573 n 



4/3^1/3^2/3 



An 



4/3 



(32) 



where A is the expression for the largest Lyapunov exponent 
describing the divergence of two initially close fluid particles, 
Th is Lyapunov exp onent was earlier obtained for such flows 
by iTuritsvnl (|2007T ) . In terestingly, this is exactly the s ame 
scaling with shear as in lVishniac fc Brandenburg i|1997l ). In 
the absence of shear the small-scale dynamo can still operate 
llChertkov, Kolokolov and Vergassolalll997l ') with j n ~ n 2 . 

IProctorl (120071 ) have also considered a model similar to 
ours, although somewhat simpler and more relevant to the 
solar dynamo, using multiscale expansions. After averaging 
over the fluctuating a effect, he still finds an effective a effect 
from which he obtains a dynamo which grows in the mean 
(as opposed to mean-square in our case). His results give 
the scaling, fc p eak ~ S and ■y ^ S 2 which disagree with the 
DNS results of lYousef et all IJ2008h . 
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APPENDIX A: DERIVATION OF EQUATIONS 

(7) AND (8) 

It is possible to derive ([8| via a method described by 
iBrissaud fc Frischl (|l974T ). This method is superior to the 
one described in Appendix iBl in the sense that this can be 
applied even when a is not necessarily white in time, but 
has (small) non-zero correlation time. On the other hand 
it is more cumbersome to apply this method to calculate 
the higher order moments of the PDF of the magnetic field. 
For the sake of completene s s we r eproduce below the calcula- 
tions of lBrissaud fc Frischl II 19741 ) as applied to our problem. 
Let us write symbolically the evolution equations for 
the mean-field in Fourier space in the following way, 



d t B = [M d + M] B. 



(Al) 



Here, B = (B x ,B y ) is a column vector, M d is the deter- 
ministic part of the evolution equations (i.e., the part that 
depends on rjij), 



k f) yv 

O —K Tjyy 



(A2) 
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which is also independent of time, and M is the random part 
(i.e., the part that depends on ohj), with 



M = ik 



O^xx 



ikA. 



(A3) 



To begin, we do not assume that a is white-in-time but that 
it has finite correlation time T CO rr. Later we shall take the 
limit of T corr — >• in a specific way to reach the white-noise- 
limit. This is equivalent to the regularization in Appendix 

m 

In this section, for simplicity, let us choose our units 
such that at t = 0, B — (1, 1). In that case the solution to 
(|A1|) can be easily recast in the integral form 



B(t) = e + 



di'e^'-^M (*')#(*') 



(A4) 



Note that we have (M) = because we have assumed all the 
components of a to have zero mean. Iterating this equation, 
we obtain 



Bit) 



+ I dt' 



j.l M d (t-t'). ./,/n M d t' 

at e k 'M(t )e 



dt"e Md{t - t ' ) M{t')e Md{t '- t " ) M{t")B{t" 



(A5) 



+higher order terms. Let us also assume that the strength 
of the fluctuations of a are finite and bounded by a. Equa- 
tion ()A5|) is then an expansion in powers of at. To obtain 
closed equations for the first moment of the magnetic field, 
we average (|A5|) over the statistics of a and then take the 
derivative with respect to t. Remembering our earlier nota- 
tion C 1 = (B), we obtain 



d t C 1 {t) = M d C 1 (t)+ / dt'> 
Jo 



M(£)e 



M d (t-t') 



M(t')B{t')).{A6) 



This obviously is not yet a closed equation. To obtain a 
closure, note that for 



o-(t'-s) < 1, 
from (|A4[) . we have 

Bit') « e Md(t '- s) B(s) + 0(a(t' - s)). 



(A7) 



(A8) 

Substituting (|A8|) in the integrand of the double integral in 
161) we obtain the factorization, 



(A9) 



(M(i)e Md(t - t ' ) M(£')e Md(t '- s) e( S ))^ 

<M(t)e Md(t -'' ) M(f'))<e Md(t '~ s) S( S )), 

if we assume that 

t'-s>T corr . (A10) 

Equations (|A7|) and (|A10|) can both hold only for small Kubo 
number, 



K = (7 T cor r "C 1. 

Substituting (|A9|) in (|A6|) and using again 



(All) 
we obtain 



d t C x {t) = M d C 1 (t)+ / dt'{W(t)e ud{t - t ' ) W{t'))C 1 (t'), (A12) 



Following iBrissaud fc Frischl (|1974T ). we shall call this equa- 
tion the Bourret equation. For small Kubo number we can 
invert (IA8I) to have 



B(t') 



-M d (t-t') 



B(t). 



Averaging (|A13|) over the noise we obtain 
C\t') « e-^^-^C^t). 



(A13) 



(A14) 



Substitute this back into (|A12|) . noting in addition that for 
short-correlated a and t 3> T corr , we can replace r = t — t' in 
(|A12|) and extend the integral from zero to infinity to obtain 



dtC 1 = M d C 1 + 



dr(M(r)e M T M(0)e~ K ' 



)C\ (A15) 



where we have omitted the t argument on all C x (t) for 
brevity. To get this result, remember that the integral above 
gives negligible contribution for r 3> Tcorr and the matri- 
ces M d and M(r) do not necessarily commute. To go to 
the white-in-time limit we need to take the limit a — > oo, 
Tcorr — > in such a way that the product o- 2 T CO rr remains 
finite. In this limit the Kubo number goes to zero and the 
various approximations made above become exact. In this 
limit the integral in (|A15P reduces to 



dtC 1 ^) = M d C 1 (i) + (MM)C 1 (i) 
= [M d - fc 2 (AA)] C x (t). 



(A16) 



The correlator on the right hand side of (|A16f) can be ob- 
tained by using (|A3f) together with (|3}- This reproduces J§J. 
If instead of @, (|20p is used, (|2ip can be obtained. 

Finally, note that for an a effect which has finite-time- 
correlations (instead of white-in-time) higher order terms in 
the expansion in (|A5|) are needed. In such cases it may not 
even be possible to obtain closed equations like (0. 



APPENDIX B: AVERAGING OVER GAUSSIAN 
NOISE 

We explain here another technique used to derive ([TJl 
and (|8} . Let us begin by considering Gaussian vector- valued 
noise Vj{t) (not necessarily white-in-time) and an arbitrary 
functional of that, F(v). Then, 



(fM"i(t))= dt'{u,(t)Mt')) 



8F 
5v k (t>) 



(Bl) 



where the average factorises by virtue of the Gaussian prop- 
erty of the noise. Here the operator 5{-)/5vk is the func- 
tional derivative with respect to v. This useful identity often 
goes by the name Gaussian integration by parts; see, e.g., 
Zinn- Justin! (I1999 | ). Section 4. 2 for a pro of; see also e.g., 
Frischl (Il996l ). iFrisch fc Wirthl (|l997n . or iMitra fc Pandfil 
l|2004h . where this method has been used to derive closed 
moment e quations for the K raichnan model of passive scalar 
advection (|Kraichnanl 1 19681 ) . 

To obtain an evolution equation for C = ((B x ), (B v )) 
we average each term of (0 and (O over the statistics of 
Qij . Terms which are product of components of ay and B x 
or By can be evaluated by using the identity in 
particular, 



Bill. In 



\(%yx-Dx ) 



D h 



*W(*)<MO>(^ 

SB x {t) 
x 6ctki(t)' 



(B2) 
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Here we have considered the magnetic field to be a func- 
tional of otij. Substituting the covariance of etij from @, 
integrating the 5 function over time and contracting over 
the Kronecker deltas, we obtain the last equality in (|B2|) . 

The functional derivatives of the components of the 
magnetic field with respect to otij can be obtained by first 
formally integrating ([5]) and ((6j to obtain B x (t) and B y (t), 
respectively, and then calculating their functional deriva- 
tives with respect to a^-. We actually need the functional 
derivative 5(B x (t))/Saki(t') and then take the limit i — > i'. 
This is a non-trivial step due to the singular nature of 
the correlation function of {otki(t)oikl(t')) as £ — > t'. To 
get around the difficulty it is possible to replace the Dirac 
delta function in ([3| with a regularise d even funct i on an d 
then take limits. We refer the reader to IZinn- Justin! |l999), 
Section 4.2, for a detailed discussion. This regularization is 
equivalent to using the Stratanovich pr escription for th e set 
of coupled SDEs © and ©; see, e.g.. lGardinerl (|l994T ). For 
reference, all the non-zero functional derivatives needed are 
given below: 



SB x (t) _ 
5a yx (t) 
SB y (t) 



-ikB x , 



= ikB x 



5B x {t) 
Sa vv (t) 



~ikB y , 



P^=ikB y . 

5a X y(t) 



8a xx (t) -"'■ < - -"■ (B3) 

In particular, since there is no a yy term in ((BJ, we have 

S(B y (t))/Sa yy (t) = 0,so 

SB V 



(OlyyDyj — L)y 



Sa v 



0. 



(B4) 



Putting everything together we can now average (J5J) over the 
statistics of a to obtain, 

dt{B x ) — — (ika yx B x ) — (ika yy B y ) 

+rj yx k 2 {B y ) ~rj yy k 2 {B x ) 
— —ik(—ik)D yx {B x ) 

+rj yx k 2 (B y ) - Vyy k 2 (B x ) 
= ~(%w + D yx )k 2 (B x ) +riy X k 2 (B y ). (B5) 

This gives us the first row of the matrix Ni in (|SJ). Applying 

the same technique to ((6j, the second row can be obtained. 

The same technique can be applied to obtain the matrix 

N2. Here we end up with evaluating terms of the general form 



(a.ijB k Bi) — Di 



5B k ~ t ~ SB? 
k Bi+B h 



5a 



' Soiii 



(B6) 



The only non-zero functional derivatives are given in (1B3I) . 
We show the calculation explicitly for the two following ex- 
amples: 



(a yx B x B x ) = D y 



SB X -„ SB* A 
B* X + ^^B, 



D yx [~ik(B x B* x ) + ik{B x B*)] = 0, (B7) 



and 

\ a yyB x B y ) = L) yy 



SB y _ §: + 6Bl A 



So 



8a v 



ikDyy(ByB*y). 



(B8) 



Note further that, instead of using the functional cal- 
culus above, the same evolution equations for C 1 can be 



obtained by using the technique due to iBrissaud fc Frischl 
|l974J). This is demonstrated in Appendix 1X1 



APPENDIX C: A ZERO-DIMENSIONAL 
MEAN-FIELD MODEL WITH FLUCTUATING a 

A simpler mean-field model in a one-mode trunca- 
tion, but with fluctuating a effect, was introduced by 
IVishniac fc Brandenburg l| 19971 ): see J25} and lf26)l. For this 
model we define the following moments of successive orders, 

C 1 = ((B X ),(By)), (CI) 

C 2 = ((Bl),(B 2 y ),(B x By)), (C2) 

C 3 = {(Bl),(B 2 x By),(B x B 2 v ),(B 3 y )), (C3) 

C 4 = ({B 4 x ),{B 3 x By),{B 2 B 2 y),{B x B 3 y ),{B y )). (C4) 

Each of these moments satisfies a closed equation of the form 

d t C p = N P C P . (C5) 

The matrices N p can be found by applying the technique 
described in Appendix [B] and by using the covariance of a 
given in (|27[) . We give below the first four matrices: 



Ni = 



N 3 = 



N 4 = 



>h 











S 


-Vt 


5 






2 Vt 


2D 


" 









-2r, 


t -25 


: 




-s 





-2*7t . 






-3??t 





6D 


" 




-S 


-377 


t 


D 







-2S -3% 





) 








-3S 


-3r?t 




-4rj t 





12D 








-s 


—477 


t 


6D 








-2S -4?jt 





2D 








-3S 


-4jjt 














-45 


—4?7 



(06) 



(C7) 



(C8) 



(C9) 



The growth rate at order p is defined to be C p ~ exp(p7 P i). 
This gives 71 = — rjt, i.e., there is no dynamo. But this also 
gives dynamo modes with positive eigenvalues given by 

lp ~S 2/3 Dli 3 , p = 2,3,... (CIO) 

The s ame result was obtained bv lVishniac fc Brandenburg! 
I|1997h for 72 by using a different method. 

Note the striking similarity between matrix N2 in (IC7|I 
and matrix N2 in (|10[) . A trivial way of generalising (IC7|I 
to one spatial dimension is to replace D and rj t in (IC7|I 



by k 2 D and fc 2 ?7t, respectively. The solution of the resultant 
eigenvalue problem gives the scaling, 7 ~ S and fc pcak ~ </S, 
Thus, (|C7|I for this zero dimensional model is equivalent to 
(|10|l in the space-time model. 
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